INTEGRABLE KP AND KDV CELLULAR AUTOMATA 
OUT OF A HYPERELLIPTIC CURVE 



MARIUSZ BIALECKI 



Abstract. The goal of this paper is to present a solution of the cellular 
automaton associated with the discrete KdV equation using algebro-geometric 
solution of the discrete KP equation over a finite field out of a hyperelliptic 
curve. 



1. Introduction 

Cellular automata (CA) are (solutions of) completely discrete dynamical sys- 
tems, for which values of both independent and dependent variables are discrete. 
The concept of CA emerged in interaction with computational machines and hence 
CA are convenient tools for computer simulations of various phenomena. There is 
also an approach to CA, in which the main objective is to obtain "analitic" solu- 
tions (or full solutions or in a form good for investigation of their global properties) 
without need of performing step by step calculations. This leads naturally to notion 
of integrability and integrable cellular automata (ICA). (For references to ICA see 

ao 

Recently a new method of construction of ICA was proposed in Its main 
idea is to keep the form of a given integrable discrete system and to transfer the 
algebro-geometric method of construction of its solutions ^01 Q] from the complex 
field C to a finite field case. In this framework there were constructed finite field 
versions of multisoliton solutions for the fully discrete 2D Toda system (the Hirota 
equation) in [gj and for discrete KP and KdV equations (in Hirota form) in [4] and 
also algebro-geometric solutions of the discrete KP equation out of a hyperelliptic 
curve in [H]. 

In this paper we extend previous works to construct algebro-geometric solution 
of the discrete KdV equation out of a hyperelliptic curve. To reach this aim we 
obtain a solution of dKP equation in a form compatibile with the reduction from 
the dKP equation to the dKdV equation. 

The paper is constructed as follows. In Section [21 we first summarize the finite 
field version of the algebro-geometric construction of solutions of the discrete KP 
and KdV equations. In Section^lwe apply the method to construct solution of the 
discrete KP and KdV equations starting from an algebraic curve of genus two. 
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2. The finite field solution of the discrete KP equation out of 
nonsingular algebraic curves 

We first shortly recall algebro-geometric construction of solutions of the discrete 
KP equation over finite fields and its reduction to the discrete KdV equation j^lQlini- 
An algebro-geometric approach in case of complex field C is described in ^J. An 
algebraic curves over finite fields are exposed for example in [151 [T7| . 

2.1. General construction for the dKP equation. For the general construc- 
tion we need an algebraic projective curve C/K (or simply C), absolutely irreducible, 
nonsingular, of genus g, defined over the finite field IK = Fq with q elements. By 
C(K) we denote the set of K-rational points of the curve. By K denote the algebraic 
closure of K, i.e., K = Ufci '^q'^^ ^ind by C{K) denote the corresponding infinite set 
of K-rational points of the curve. Denote by Div(C) the abelian group of the divisors 
on the curve C. The action of the Galois group G'(K/IK) (of automorphisms of K 
which are identity on K) extends naturally to action on C(K) and Div(C). A field of 
K-rational functions on the curve C we denote by K(C) and the vector space L{D) is 
defined as {/ G K(C) | (/) > ~D}, where D € Div(C) and (/) = Y.PeC ordp{J) ■ P 
is the divisor of the function / G K(C). 
On the curve C we choose: 

(1) four points G C(K), i = 0,1,2,3, _ 

(2) effective K-rational divisor of order g, i.e., g points G C(K), 7 = 1, . . . , 5, 
which satisfy the following K-rationality condition 

VctgG(K/K), a{B^)=B^,. 

We assume that all the points used are distinct and in general position. In particu- 
lar, the divisor Y^l^=i ^1 {and a divisor _D(ni, 712,713) defined below) is non-special. 

Definition 1. Fix K-rational local parameter to at Aq. For any integers ni, 722, 713 G 
Z let divisor D{ni, n2, n^) be of the form 

g 

D{ni,n2,n3) = ni{Ao - Ai) + n2{Ao - A2) + usiAo - A3) + ^2^^- 

7=1 

The function ^/;(ni, 7^2, ^3) (called a wave function) is a rational function on the 
curve C with the following properties 

(1) the divisor of the function satisfy (V') > —D, i.e. ip G L{D), 

(2) the first nontrivial coefficient of its expansion in to at Aq is normalized to 
one. 

Existence and uniqueness of the function ip{ni,n2,n3) is due to application of 
the Riemann-Roch theorem with general position assumption and due to normal- 
ization. Moreover, the function ip{ni,n2,n3) is K-rational, which follows from 
K-rationality conditions for sets of points in their definition. 

Remark. Notice that the function ^/'(ni, rt2, ^3) has g zeros not explicitely specified 
in the Definition 

The next step of the construction is to obtain linear equations for wave functions. 
The full form of such equation is in case when the pole of ip{ni,n2, ns) at Aq is of 
the order exactly (ni + n2 + ns) and respective zeros at Ai are of the order m, for 
i — 1,2, 3. We will call this case generic. Having fixed K-rational local parameters 
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ti at Ai, i — 1,2,3, denote by C^*^ ("i, ^3), « — 1,2,3, the K-rational coefficients 
of expansion of -(/'('t-i, ^2, ^3) at Ai, respectively, i.e.. 



t/;(ni,n2,n3) =<f'^Ci*^(ni,n2,«3)if, * = 1,2,3. 

k=0 



Denote by Ti the operator of translation in the variable n^, i — 1, 2, 3, for example 
T2ip{ni, 712, "■3) — 4'ini,n2 + 1, ^3). The full linear equation is of the form 

(1) t.tA - T,^ + = 0, z^j, = 1,2,3. 

Co 

It follows from observation, that TiTp — TjTp £ L{D), hence must be proportional to 
wave function ?/;. Coefficients of proportionality can be obtained from comparison 
(the lowest degree terms) of expansions of left and right sides of Q at the point 
A,. 

Remark. When the genericity assumption fails then the linear problem (QJ degen- 
erates to the form Tiip — ip or even to = 0. 



(2) -^-~-%r^ ^'J'^i'^^s. 



Notice that equation iQJ gives 

~ A. 
so So 

Define 

(3) p, = (_l)S.<.".C^'), * = 1,2,3, 

then equation ^ implies existence of a K- valued potential (the r-function) defined 
(up to a multiplicative constant) by formulas 

(4) — z = 1,2,3. 

T 

Finally, equations (QJ give rise to condition 

rr.\ ^2Pl T3P1 T3P2 

(5) \ = 0, 

pi pi P2 

which written in terms of the r-function gives the discrete KP equation called 
also the Hirota equation 

(6) (Tit) (T2T3T) - (r2T) {TaTir) + {T^t) (T1T2T) = 0. 

Remark. Equation 10 can be obtained also from expansion of equation iQJ at , 
where A; = 1, 2, 3, fc 7^ i, j. 

Absence of a term in the linear problem (QJ refiects, due to Remark above, in 
absence of the corresponding term in equation This implies that in the non- 
generic case, when we have not defined the r-function yet, we are forced to put it 
to zero. 
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2.2. Reduction to the dKdV equation. The discrete KdV equation [71 

(7) (Tir) r - (Tg- V) {T^T^r) + (Tgr) (TiTg'V) = 0, 
is obtained from the discrete KP equation by imposing constraint 

(8) T2Ts,T = 7T, 



where 7 is a non-zero constant. Algebro-geometric solutions of the discrete KdV 
equation can be constructed using following facts (see 0). 

Lemma 1. Assume that on the algebraic curve C there exists a rational function 
h with the following properties 

(1) the divisor of the function is (h) = A2 + — 2Aq, 

(2) the first nontrivial coefficient of its expansion in the parameter to at Aq is 
normalized to one. 

Then the wave function ip satisfies the following condition 

(9) TsTaV' = hi:. 

Remark. Existence of such a function h implies that the algebraic curve C is hy- 
perelliptic. 

Proposition 2. Let h be the function as in Lemma^ Assume additionally that 

(10) h{Ai) = 1. 

Denote by 82 and 5^ the respective first coefficients of local expansion of h in pa- 
rameters t2 and tz at A2 and A3, i.e. h = t2{S2 + ■ . .), h = t3{S3 + . . .). Then 
the function 

(11) f ^ ^5-"=("^-l)/2(_J3)-"3(«3-l)/2 

satisfies the discrete KdV equation ijjjl. 

3. A "hyperelliptic" solution of the discrete KP AND KdV equation 

Our goal here is to construct a solution of the dKP equation to which we can 
apply the reduction scheme described above. We are forced to perform steps of the 
construction (see also [2] for details) starting from a hyperelliptic curve. In detail 
we deal with a curve of genus g = 2 but the technical tools used here can be applied 
directly to hyperelliptic curves of arbitrary genus. 

3.1. Hyperelliptic curves and Jacobian picture of the construction. In the 

following we use an afRne picture of a hyperelliptic curve. It is motivated by the 
fact that general hyperelliptic curve can be tranformed to the form with only one 
point at infinity (see ^^). 

Definition 2. A hyperelliptic curve C of the genus g over a field K is given by 

(12) C: v"^ + h{u)v- f{u) = 0, 

where h{u) G K.[u] is a polynomialof degree at most g, f{u) G K[u] is a monic 
polynomial of degree 2(7+1, if there is no points {u = x,v = y) e K x K which 
satisfy equation ltT2ll and equations 2y + h{x) — and h'{x)y — f'{x) = 0. 
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By P we denote the point opposite to P, i.e. conjugate with respect to hyper- 
elliptic isomorphism. Denote by Div°(C;K) the abelian group of the K-rational 
divisors on the curve C and by J{C; K) the group of eqivalence classes of K-rational 
degree zero divisors Div°(C;IK) modulo the K-rational principal divisors, i.e. divi- 
sors of a functions K(C). (In terms of algebraic geometry J(C; K) is identified with 
the group of K-rational points of the Jacobian of the curve C |12l 114) . 1 

Two divisors A,B ^ Div"(C; K) are equivalent (we write A ~ B) if B = A + (/) 
for some function / S K(C). The class of a divisor A in the divisor class group 
J(C;K) is denoted by [A\. For a hyperelliptic curve C (of genus g) each class 
of equivalence [A\ G J(C;K) has a unique representant of the form of reduced 
divisor 13j. 

Definition 3. A divisor D e Div°(C;K) of the form 

D ^^X^ ~k ■ Ao, 

7=1 

where G C \ {Aq} is called reduced if 

(1) ^<5 

(2) X^ ^ X^. for all 7 7^ 7' 

Let us present in this picture the description of the wave function ■0 and of the 
T-function. Consider the following divisor D{ni,n2, n^) e Div'^(C; K) of degree zero 

a 

D{ni,n2,n3) = ni{Ao - Ai) + ^2(^0 - A2) + 713(^0 - A3) + - .9 • Ao, 

7=1 

with linear dependence on ni, n2 and n^. Its equivalence class in J(C;K) has the 
unique K-rational representant of the form of reduced divisor 

fe 

X{ni,n2,n3) = ^ X^(ni, n2, 713) - k ■ Aq. 
7=1 

This equivalence is given by a function whose divisor is 

k g 

ni[Ai - Ao) + 722(^2 - Ao) + 713(^3 - Ao) + ^ X^{ni,n2,m) -^B^ + {g- k)Ao. 

7=1 7=1 

If we normalize such a function at Aq acording to Definition it becomes the wave 
function ip. Notice that if fc ^ g, it means that the pole of the wave function 
at ^0 is of the order less then {ni + n2 + ns), and it is a non-generic case, thus 
T(ni, 712,713) = 0. 

Remark. Points Xj, indicate zeros of the wave function which are not explicitely 
specified in the previous construction. 

3.2. A curve and its Jacobian. Consider a hyperelliptic curve C of genus 5 = 2, 
defined over the field F7 and given by the equation 

(13) C : + uv = u^ + 5u'^ + 6u^ + U + 3. 

The {u,v) coordinates of its Fy-rational points are presented in Tabled 

The only two special points of the curve are (6,4) and the infinity point 00. 
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oo 


Po 


1 


(1,1) 


(1,5) 


2 


(2,2) 


(2,3) 


3 


(5,3) 


(5,6) 


4 


(6,4) 


P4 



Table 1. Fy-rational points of the curve C. The point opposite 
to P (conjugate with respect to the hypereUiptic automorphism) 
is denoted by P. 



i 




P^ 


pr 


Pr 


5 


(0,21) 


(0,28) 


P5 


P5 


6 


(3,9) 


(3,44) 


Pe 


Pe 


7 


(4,26) 


(4,33) 


P7 


Pj 


8 


(7,5) 


(7,44) 


(42,5) 


(42,9) 


9 


(8,22) 


(8,26) 


(43,29) 


(43, 33) 


10 


(11,5) 


(11,47) 


(46,5) 


(46, 12) 


11 


(12,6) 


(12,45) 


(47,6) 


(47, 10) 


12 


(13,14) 


(13,29) 


(48,35) 


(48, 22) 


13 


(14,8) 


(14,34) 


(35,43) 


(35,27) 


14 


(15,13) 


(15,28) 


(36,48) 


(36,21) 


15 


(16,17) 


(16,23) 


(37,38) 


(37, 30) 


16 


(17,0) 


(17,39) 


(38,0) 


(38,18) 


17 


(18,4) 


(18,41) 


(39,4) 


(39,20) 


18 


(19,9) 


(19,28) 


(40,44) 


(40,21) 


19 


(20, 12) 


(20,31) 


(41,47) 


(41,24) 


20 


(22,4) 


(22,30) 


(29,4) 


(29,23) 


21 


(25,6) 


(25,32) 


(32,6) 


(32,25) 


22 


(27,7) 


(27, 22) 


(34,42) 


(34, 29) 



Table 2. F49-rational points of the curve C (which are not F7- 
rational); P'^ denotes conjugate to P with respect to the action of 
the Frobenius automorphism. 



We identify the field F4g with the extension of F7 by the polynomial + 2, 
i.e., F49 = F7[a;]/(a;^ + 2). It is convenient to introduce the following notation: the 
element k G F49 represented by the polynomial /3a; + a is denoted by the natural 
number 7/3 + a. The Galois group G'(F49/F7) — {id, a}, where tr is the Frobenius 
automorphism, acts on elements of F49 \ F7 in the following way: fc = 7/3 + a 1— > 
<7{k) = 7(7 — (3) +a. The coordinates of F49-rational points of the curve (which are 
not F7-rational) are presented in Table |2| 

The full description of the group J(C;F7) is given in Table EI where we have 
choosen as a reference point Aq the infinity point 00. The divisor Di = Pi — 00 
generates the subgroup of order 31 and the divisor D4 — P4 — 00 generates the 
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n 


[nDi]r 


9oin) 


[nDi + 1)4] r 


.91 ("-) 


n 
U 







1 


(6,4) 


— 00 


-1 
1 


1 


(1,1)- 


oo 


1 


(1,1) + (6,4) 


- 2oo 


5+5u+3u^+v 


6+4u+ti^ 


2 


(1,1) + (1,1)- 


2oo 




(12, 45) + (47, 10) 


- 2oo 


l+5u^+v 


{2+uy 


2+5u+u'^ 


3 


(5, 6) + (5,6) - 


2oo 


l+u+4u'^+v 
(2+u)(5+«) 


(15,28) + (36,21) 


- 2oo 


6u^ -i-v 

2+M^ 


4 


(2, 3) + (5,3)- 


2oo 


2+Au'^+v 


(7,44) + (42,9) 


- 2oo 




5+4m+m^ 
iu+2u^+v 


4+61/ + 1/- 


5 


(19,9) + (40,44) - 


2oo 


(11,5) + (46,5) 


- 2oo 


6+()n-n- + u 




3+6u+u^ 


6 


(22, 4) + (29,4) - 


2oo 


5+2«+6ti^ + t! 


(18,41) + (39,20) 


- 2oo 


5+3u+5u^+v 


(2+tt)(5+«) 
5+6u+2u^ +v 


5+3u+u-^ 


7 


(2, 3) + (5,6)- 


2oo 


(16, 17) + (37, 38) 


- 2oo 


5+4«+4u^+ii 


5+2«+«^ 




8 


(27, 22) + (34,29) - 


2oo 


l+3M+2tt^+t> 


(17, 39) + (38, 18) 


- 2oo 


3+2ii+«^+t) 






(5+M)(6+ti) 


9 


(14, 34) + (35,27) - 


2oo 


(l+u)(5+ti) 
3+5u+5u^+v 

(5+t/.)' 


(1,5) + (2, 2) 


- 2oo 


6 + u 


10 


(2, 2) + (6,4) - 


2oo 


(2,2) ■ 


— 00 


1 


11 


(2, 3) + (2,3) - 


2oo 


3+2m+m^ 


(1,1) + (2, 2) 


- 2oo 


4+2m^+i> 
3+5M+ti^ 


12 


(13, 14) + (48,35) - 


2oo 


3+6u+4u^+v 


(8, 22) + (43, 29) 


- 2oo 


2+4u+v 


2+2m+«^ 


(2+m)(6+«) 


13 


(20, 12) + (41,47) - 


2oo 


(l+u)(2+u) 

6+5M+2tt^+t> 


(1,5) + (5,3) 


- 2oo 


6 + u 


14 


(5, 3) + (6,4)- 


2oo 


(5,3) 


— 00 


1 


15 


(25,32) + (32,25) - 


2oo 




(1,1) + (5, 3) 


- 2oo 






(2+ti)(6+«) 


16 


(25, 6) + (32,6) - 


2oo 


6+5u+2u^+v 

(l+tt)(2+ii) 


(1,5) + (5, 6) 


- 2oo 


+ u 


17 


(5,6) + (6,4) - 


2oo 


5u-\-u'^ -\-v 
2+2u+u'-' 


(5,6) 


— 00 


1 


18 


(20, 31) + (41,24) - 


2oo 


3+2m+«^ 


(1,1) + (5, 6) 


- 2oo 


2+4u+t) 
3+5M+ti^ 


19 


(13, 29) + (48,22) - 


2oo 


{5+uf 


(8, 26) + (43, 33) 


- 2oo 


4+2M^+t> 
(5+u)(6+«) 


20 


(2, 2) + (2,2) - 


2oo 


^ '2 


(1,5) + (2, 3) 


- 2oo 


6 + U 


(l+ti)(5+«) 


21 


(2, 3) + (6,4) - 


2oo 




(2,3) 


— 00 


1 


22 


(14, 8) + (35,43) - 


2oo 


1+3u+2m^+i; 

5+2«+«2 


(1,1) + (2,3) 


- 2oo 


3+2u+u^+'!) 
3+it+u2 


23 


(27, 7) + (34,42) - 


2oo 


5+6M+2ti^+r 


(17,0) + (38,0) 


- 2oo 


5+4M+4ti^+ii 


(2+«)(5+«) 


5+3ti+u^ 


24 


(2, 2) + (5,3)- 


2oo 


5+2w+6m^+i; 


(16, 23) + (37, 30) 


- 2oo 


5+3m+5m^+i; 




3+6u+u^ 


25 


(22,30) + (29,23) - 


2oo 


4u+2u'^+v 


(18, 4) + (39, 4) 


- 2oo 




5+iu+u'^ 


4+6u+u2 


26 


(19, 28) + (40,21) - 


2oo 


2+4«2+D 


(11, 47) + (46, 12) 


- 2cx) 




(2+u)(5+u) 


2+1i2 


27 


(2, 2) + (5,6) - 


2oo 


l+tt+4M^+t> 


(7, 5) + (42, 5) 


- 2oo 




(2+«)- 


2+5u+u2 


28 


(5, 3) + (5,3) - 


2oo 


tt+5ti 
(6+tt)^ 


(15,13) + (36,48) 


- 2oo 


6+4ti+tt^ 


29 


(1,5) + (1,5)- 


2oo 


(6 + w) 


(12, 6) + (47,6) ■ 


- 2oo 


5+5«+3u^+r 


(1+«)(6+m) 


30 


(1,5)- 


oo 


(6 + w) 


(1,5) + (6,4) 


- 2oo 


6 + u 



Table 3. The group J(C; F7) as the simple sum of its cydic sub- 
groups; [X]r denotes the reduced representsnt of the equivalence 
class of the divisor X; gQ[n), gi{n) — transition functions (see 
main text). 
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subgroup of order 2. For n e {0,1,..., 30} and m e {0, 1} we present the reduced 
representants of elements [nDi + mZ?4],. of J{C;¥j) and the transition functions 
gm{n) which are given by the equation 

[nDi + mD4]r + Di = {g„i{n)) + [{n + l)Di + mDi]r, 

and normaUzed (numerators and denominators are monic polynomials). We will 
use them in the construction below. 

3.3. Construction of the wave and t functions. In order to find a solution of 
the discrete KdV equation let us fix the following points of the curve C: 

Ao^oo, yli = (2,2), A2 = (l,5), A3 = (1,1), 

with the uniformizing parameters io = u^/w, ti = u — 2, t2 — — u — I, and 

Si = (12,6), 52 = (47,6). 

Then 

Ai - Ao ^ lODi + Di, A2 - Ao -D^ ^ 30Di, A3 - - A, 
B1 + B2- 2Ao - 29Di + Di, 

and the points Xi(ni, 712, 713) and X2(ni, 712, 713), where the wave function V'("ij "2, ^3) 
has additional zeros (here Xi can be 00) can be found from Table Eland 

(14) Xi{ni,n2,n3) + X2('^i, "2, "3) - 2oo = [nDi +mD4]r, 
where 7t e {0, 1, . . . , 30} and m £ {0, 1} are given by 

(15) n EE 29 - (713 - 712 + lOni) mod 31, 

(16) m = 1 — ni mod 2. 

Remark. The choice of the infinity point 00 as Aq is a violation of the assumption 
of general position of points used in the construction (00 is the Weierstrass point of 
the curve C). This will not destroy the construction but in some situations, which 
we will point out, will affect uniqueness of the wave function. We remark that such 
a choice is indispensable in reduction of the method from the discrete KP equation 
to the discrete KdV equation (see, for example ^UBI)- 

The key idea in constructing of the wave function is to express 7/1(711,772,713) 
for any parameters from (771,772,773) G by a set of functions related to J(C;IK) 
(transition functions and few auxiliary functions). Let us introduce functions hi 
and /i4 corresponding to generators of the two cyclic subgroups of J{C;¥j). The 
function hi with the divisor 31I?i ^ and normaHzed at the infinity point is equal 
to 

30 

hi = ]J.9o(i): 

1=0 

and reads 

/7l = 1 + 2li + + 47t^ + 37t^ + 3u'^ + 47t^ + 

4w^" + 2u" + 5u^^ + 2u^^ + ^y}^ + 37t^^ + (Sii + 2v?+ 

hu^ + 4m''^ + 677® + 47i^ + 37i^ + 57il" + 57i" + 47*12 _^ yl3\^ 
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where we also used equation of the curve l(T3ll to reduce higher order terms in v. 
The normaUzed function /14 with the divisor 2D4 is 

= u — 6. 

Let us introduce other auxilUary function /i and /2 to factorise the zeros at Ai and 
A2 of the wave function. Notice that 

(2, 2) + 21(1,1) + (6, 4) -23(30-0, 

which impUes that there exists a polynomial function on C with simple zero at Ai 
and other zeros in the distinguished (by our choice of description of J(C; F7)) points 
(1, 1) and (6,4). Define /i as the unique such function normalized at the infinity 
point 00, then 

fi = l + 5u + u^ + 4:u^ + 6u^ + 4u^ + 4w^ + + 4u^+ 

6m" + (6 + 4u + 2^2 + 5u^ + 6u^ + 6u^ + + + u^)v. 
The zeros at A2 can be factorised using function 

/2 - {U - 1). 

Uniqueness of the wave function ijj implies that it can be decomposed as follows 

(17) ii){ni,n2, ns) = ^p^q W(mi,m2), 
where new variables mi i ni2 are given by 

(18) 21ni — ?i2 + rt3 = 'Hp — mi, toi S {0, 1, . . . , 30}, 

(19) ni = 2(7- 7712 , 7772 6 {0,1}, 

and the function W{mi,m2) has the divisor 

(20) miDl + 7772134 + ^1 (™l , 7772) + F2 (mi , 7772) - (12, 6) - (47, 6). 

The additional zeros 

11(7771,7772) +l2(mi,?772) = ^1(771,772,773) + ^2(771, 772, 773), 

can be found by projection from ?7-variables into the 777- variables. 

To find the functions W{mi,m2) for all ttii G {0, 1, . . . , 30} and 7772 G {0, 1} let 
us notice that W{0, 0) = 1 and W{0, 1) is given by 

iW{0, 1)) = Di + (1, 5) + (1, 5) - (12, 6) - (47, 6), 

and hence can be written in a form 

2 + 377 + 4.^ + . ^ 

Define the multipliers Wm2(j'^i) as follows 

iy(777i, 7772) = Wm2i'm'i)W{mi - 1, 7772), 
for 7771 g {0, 1, ... , 30} and 7772 £ {0, 1}. Equations (1111-1119 and gives 

Wm2(r77i) = gm{n), 

where 

7772 = 1 — 777 mod 2, 777i = 29 — 77 mod 31. 
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5 10 15 20 25 30 35 5 10 15 20 25 30 35 



■ - (00), ■ - (01), -(02), -(03), ■-(04), -(05), -(06). 

Figure 1. The Fy-valued solution of the discrete KP equation 
out of genus 5 = 2 hyperelliptic curve C; n2 = —1, 0, 1, 2 for sub- 
seqent figures, ni = 0, 1, . . . , 34 (horizontal axis), — 0,1, ... , 34 
(vertical axis). 



Setting 10^2(0) = W{0,m2) we obtain 

mi 

W{mi,m2) = J|w,„2(i). 

1=0 

Together with factorisation ifTTjl it gives the wave function ip for all (ni, ^2, "3) £ 

Remark. For (mi, 7712) = (29, 1) we have Xi = X2 = 00. Because the infinity point 
00 is the Weierstrass point of order two, there exist functions with divisor of poles 
equal to 2oo. This means that tp is not uniqely determined in this case. However it 
is natural to keep the divisor of -0, and therefore V' itself, exactly like it is given from 
the flow on J(C;K). Notice that because for Xi = X2 = 00 we stay in non-generic 
case, then this ambiguity does not affect construction of the r-function. 



INTEGRABLE KP AND KDV CELLULAR AUTOMATA 
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20 40 60 80 



■ - (00), ■ - (01), -(02), ■-(03), ■-(04), -(05), -(06). 

Figure 2. The Fr-valued solution of the discrete KdV equation 
out of genus g — 2 hyperelliptic curve C; ni = 0, 1, . . . , 79 (hori- 
zontal axis), 713 = 0, 1, . . . , 79 (vertical axis), (^2 = 0). 



The coefficients Co ('^i>"2,?^3), k — 1,2,3, of expansion of the wave function 
can be obtained from factorisation (tTzIl and are given by 

(21) (2,2): C^o^\n,,n2,ns) = 6"'5HPW{m,,m2)\t,=o, 

(22) (1,5): Co'^(ni,n2,n3) = 2"MW(mi, m2)|t,=o, 

(23) (1,1): Co ("-i,n2,n3) = 6^'4« -^^ |t3=o- 

'3 

Using definition of the r function for nonzero pi, i.e. equation and putting 
T = for points related with nongeneric case we obtain a solution of the discrete 
KP equation © taking value in the finite field F7. This r-function is presented in 
Figure 

To obtain a f -function being a solution of the discrete KdV equation we use the 
formula llll|l . For our settings we have 62 = S3 = 1, so finally 

f — ^(_l)-"3(n3-l)/2^ 

The f -function for dKdV is presented in Figure [21 
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